Non-Markovian dissipative dynamics of two coupled qubits in independent reservoirs: 
a comparison between exact solutions and master equation approaches 
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The reduced dynamics of two interacting qubits coupled to two independent bosonic baths is inves- 
tigated. The one-excitation dynamics is derived and compared with that based on the resolution of 
appropriate non-Markovian master equations. The Nakajima-Zwanzig and the time-convolutionless 
projection operator techniques are exploited to provide a description of the non-Markovian features 
of the dynamics of the two-qubits system. The validity of such approximate methods and their 
range of validity in correspondence to different choices of the parameters describing the system are 
brought to light. 
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I. INTRODUCTION 

Despite its simplicity a two-state system is of great 
significance being it exploitable to effectively describe 
many real situations. The theoretical analysis as well as 
the practical implementation of interacting or not two- 
level systems thus represents a central topic in several 
branches of modern physics ranging from high energy to 
nuclear and condensed matter physics [l|, 0| ■ During the 
last decade the interest towards two-level systems has 
further been stimulated by the fact that a qubit repre- 
sents the basic element in the context of the new ap- 
plicative area of quantum information and communica- 
tion. It has been for example shown that a fundamental 
quantum gate like the C-NOT gate can be implemented 
using dipole-dipole interacting quantum dots modeled 
as two qubits [3[. Moreover due to the development of 
new technologies, today there are several possible routes 
to the creation of what might be termed quantum bit, 
each based on a different physical system. These include 
quantum optics, microscopic quantum objects (electrons, 
ions, atoms) in traps, quan tum dots and quantum circuits 
[1 S, H 0, HI, S EE El M El ■ In describing real systems 
however it is mandatory to take into account the effects 
stemming from the presence of the surroundings. Thus 
the dissipative dynamics of two-level systems has been 
the subject of numerous paper s appea red in literature 
in the last decades [3, El El, El 111, EI H3 . Gener- 
ally speaking the research has been developed assuming 
a Markovian environment [LH [l5j . But memory effects 
are in general present and could affect quantitatively and 
qualitatively the dynamics of the small system. Unfor- 
tunately, there are no fully systematic investigations of 
non-Markovian environments. Projection operator tech- 
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niques, such as the time-convolutionless (TCL) [21( and 
the Nakajima-Zwanzig (NZ) 0, [23[ approaches, are in 
general exploited in order to perform a description of the 
non-Markovian features of the dynamics of open systems. 
On the one hand, the NZ provides a generalized master 
equation in which the time derivative of the density oper- 
ator is connected to the past history of the state through 
the convolution of the density operator and an appropri- 
ate integral kernel. On the other hand the TCL approach 
provides a generalized master equation which is local in 
time. Intuitively, one might argue that the NZ should 
work better than the TCL approach in describing the 
memory effect, since it explicitly takes into account the 
past history of the open system. Anyway there are exam- 
ples in which the exact dynamics of the open system can 
be described by means of a master equation which is local 
in time, as in the well known case of the Hu-Paz-Zhang 
generalized master equation for the non-Markovian the- 
ory of quantum Brownian motion [24], [25[ . 

In general it is not easy to establish whether one 
method is better than the other one. In fact, the perfor- 
mance of these perturbation schemes strongly depends 
on the details of the system under investigation. In this 
paper we consider two interacting two-level systems, each 
coupled to its own bosonic bath, exactly solving their dy- 
namics in a one excitation subspace. The knowledge of 
the exact dynamics is exploited to test perturbative ap- 
proaches based on TCL and NZ techniques. We show 
that, counterintuitively, the TCL approach works better 
than the NZ one, since the latter approach does not guar- 
antee the positivity of the density matrix when the cor- 
relations inside the reservoir become moderately strong. 
On the contrary the TCL approach describes all the qual- 
itative features of non-Markovian dynamics for a wider 
range of values of reservoir memory time. 

The paper is organized as follows. In Sec. II we derive 
the exact equations governing the evolution of the two- 
qubit system coupled to two independent reservoirs, in 
the case of one initial excitation, and find their exact 
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solution. In Sec. Ill the derivation of the second order 
NZ equation is presented and the features of the second 
order TCL, derived in Ref. [26[, are recalled. An extensive 
comparison among the exact, the NZ and the TCL is 
presented in Sec. IV, while in Sec. V some conclusive 
remarks are given. 

II. EXACT DYNAMICS 

A. The model 

The physical system on which we focus our attention 
is composed by two interacting two-level systems. Each 
qubit is moreover coupled to an external environment 
modeled as a bosonic bath [l6[. Assuming h = 1, the 
Hamiltonian model describing the total system can be 
written in the following form 

H = H + H T . (1) 

Here 

fl» = ^ + ^+EE« ( 2 ) 

j=l,2 k 

is the unperturbed part containing the free Hamiltonian 
of the two qubits as well as that of the two indepen- 
dent environments. The transition frequency of the two 
two-level systems, supposed coincident for simplicity, is 
indicated by coq whereas Uz (j = 1, 2) denotes the Pauli 
operator describing the j-th subsystem. The two inde- 
pendent bosonic baths are characterized by proper fre- 
quencies , and &jr' being correspondingly the 
creation and annihilation bosonic operators. 
The interaction term 

+ E f^E^^ + ^E^e) (3) 

3=1,2 \ k k J 

includes both the direct interaction between the two 
qubits, characterized by the coupling constant Q, and 
the interaction between each qubit and its respective 

bosonic bath, with coupling constants . In eq.([3]) 

cr± = !(o4 ± toy ) are j as usual, the lowering and 
raising Pauli operators. 

It is worth underlining that the Hamiltonian model 
(U) is quite versatile in the sense that it can be success- 
fully adopted to describe many different physical systems. 
In the framework of cavity QED or circuit QED 
[E S S B Hi this model can be indeed exploited for the 
description of two atoms in spatially separated cavities as 
well as of two far enough Josephson charge, flux or phase 
qubits so that it is reasonable to assume that they in- 
teract with two different electromagnetic environments. 
Model ((T|) in addition allows to study the influence of 



spurious microwave resonators within Josephson tunnel 
junctions on the coherent dynamics of a phase qubit [27| . 
In this case the environment coupled to the spurious res- 
onator (modeled as a two state system) is a phononic 
bath. 

Very recently the Markovian dynamics stemming from 
Hamiltonian (JT|) has been analyzed [TB|. In Ref. [2^] in- 
stead the TCL approach has been exploited in order to 
investigate on the non-Markovian regime. In this paper 
we exactly solve the time-dependent Schrodinger equa- 
tion confining ourselves to the one excitation subspacc. 
At the same time we adopt the NZ technique to derive 
a non-Markovian master equation for the reduced den- 
sity operator of the two coupled qubits. Having at our 
disposal both the exact dynamics and the approximate 
master equations, a comparison may be done in order 
to test the effectiveness of the perturbative approaches. 
From now on, we work in the interaction picture defined 
by Hq in which the interaction Hamiltonian reads 

H I (t) = H ( I s) +H\ Dl) +H\ D "- ) (4) 

with 

ffW=fi(4V 2) +<T«af ) (5) 

and 

+ cr o-)^ 5 «* 6 a)t e -^o-^)t. 

k 

B. One excitation time evolution 

Let us begin by looking at the exact solution of the 
time-dependent Schrodinger equation. It is easy to verify 
that the number operator 

3=1,2 j=l,2 k 

is a constant of motion. This in particular means that, 
starting from an eigenstate of N , the system evolves re- 
maining in the subspace correspondent to the same eigen- 
value n of JV. In what follows we consider the dynamics 
of the system in the subspace with one excitation, that 
means n = 1. To this end let us suppose to prepare at 
t = the two qubits in a linear superposition of states 
with one excitation and both the baths in the vacuum 
state denoted by |0^ } ) (j = 1,2) 

\m) = (a(0)|10) + 6(0)|01)) |0i 1} 0i 2) ), (8) 
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with |a(0)| 2 + |6(0)| 2 = 1. Since [H,N] = 0, at a generic 
time instant t we may write 

|^(0) = (a(i)|10)+^)|01))|0i 1) 0i 2) ) + 

+ ioo) fE^wiiW^+E^wi ^ 1 ^). 

V k k / 

(9) 

where |lj^) denotes a state of the j — th bath (j = 1,2) 
with one excitation in the mode fc and the probability 
amplitudes a(t), b(t) and cl (j = 1,2) are solutions of 
the following system of differential equations 

f d(t) = -» (n&(f) + £ fcl 



6(t) 



^a(*)+E fe 4 2) (%re 



7 (2) p i( W0 - w < 2 ))t N 



(1)* -i(w -^ 1) )t 



(10) 



4 2 >(t) = -i6(t) ff W* e -*C«o-«r')t 



(2)* 



From eas.(|10p it is easy to verify that the amplitudes 



c[ (t) formally evolve as follows: 



.,(2) 



(*) 



(*) 



/ a(t')e-^<>-^ dt' 



(2) 



} * / fe(t')e- l( "°-^ 2)) *' dt'. 



(11) 



(12) 



Inserting these formal solutions in the equations for ait) 
and b(t) we achieve 



a(t) = -iClb(t) - i SlaWfxit-t^dt 1 
b(t) = -iQa(t) - i f* b{t')f 2 {t - t') dt' 



(13) 



where the kernel fj(t — t') is given by the correlation 
function defined as 



/i(*-o =Y,\9 { k 



{j)\2 e i{w -^'){t-t') 



(14) 



that in the continuum limit becomes 

f+OO 

f j (t-t') = duJjMe*" -"**-*'), (15) 
Jo 

Jj{oj) being the spectral density of the j — th bath. 

Making use of the Laplace transform, the system (fT3|) 
becomes 



sa(s) — a(0) = —iQb(s) — a(s)fi(s) 
s b(s) - 6(0) = -ifl a(s) - b(s)f 2 {s), 



(16) 



where a(s), b(s) and fj(s) denote the Laplace transforms 
of a(t), b{t) and fj(t — t') respectively. It is thus imme- 
diate to obtain 



a{s) 



a(0)(s + / 2 (s)) -iSlb{0) 
(s + f 1 ( s ))i s + f 2 { s )) + ^ 



(17) 



6(5) = 



b(0)(s + /i(s)) -fflq(O) 
(s + /i(3))( a + / 2 ( S )) + n 2 ' 



(18) 



Once fixed the spectral densities for both baths Ji (w) and 
J 2 {io), it is quite easy to obtain the time behavior of a(t), 
6(t) and cj^(t), simply antitransforming the amplitudes 
a(s) and b(s) given by eqs. tfTT)) and (fT5|) . The results thus 
obtained will be discussed in Sec. IV. 



III. NAKAJIMA-ZWANZIG MASTER 
EQUATION 

In this section we will apply the projection operator 
techniques in order to derive a non-Markovian master 
equation for the reduced density matrix ps(t) of the two 
qubits [25| . To this end, it is convenient to introduce a 
super-operator according to 



p^Vp = Ti B P® Pb = Ps ® Pb, 



(19) 



where ps is the density matrix of the environment. The 
super-operator V projects any state of the total system 
p onto its relevant part Vp, expressing formally the elim- 
ination of the irrelevant degrees of freedom from the full 
dynamical description of the model under scrutiny. Fol- 
lowing the NZ approach we get an integro-differcntial 
equation 



(20) 



describing the reduced dynamics of the system. Here the 
memory kernel /C(i, t') is a super-operator in the relevant 
subspace. In order to discuss the reduced dynamics we 
perform a perturbation expansions of /C(i, t 1 ) with respect 
to the strength of the interaction Hamiltonian ^25]. If we 
restrict ourselves to the second order, the two relevant 
terms of /C(t, t') can be written down as 



and 



/Ci(M') =V£(t)V 



JC 2 {t,t')=VC(t)QC(t')V, 



(21) 



(22) 



where Q = I—V and C is the Liouville super-operator de- 
fined by C{t)p{t) = -i[Hj(t),p(t)]. Starting from eq.([20| 
it is possible to demonstrate that the second order non 
Markovian master equation assumes the following form 



r-p s it) = - l [H\ s \ Ps {0)]+ 



dt' Tr B {[Hxit), [Hj{t'), psit') ® pb}}} , 

(23) 



that, contrarily to the TCL master equation, is a non- 
local evolution equation. After some manipulations it is 
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possible to recast eg. (|23|) in the more compact form 
j t p s {t) ^- l [H ( I s) lPs (0)} + C^p s (t)+ 

where the dissipators are expressed by 

C^ Ps (t) = - fdt'lH^itUH^it'lpsit')]} (25) 
Jo 

and 

CV>*>p 3 {t)= fdt 1 {g j {t-t')[a^p{t l \a ( i ) ]+ 

Jo 1 (26) 

+g*(t-t'){a^,p s (t')a^}}. 

In cq. (|2"6"l) gj(t — t') is the correlation function, that in 
correspondence to a thermal bath with T = coincides 
with fj(t — t') defined in eq. (|14[) . Once again, in order to 
solve the master equation (|24|) we can exploit the Laplace 
transform. The results we have obtained are reported in 
the next section where they are compared with the exact 
dynamics as well as with the results obtained exploiting 
the TCL approach [26[. The master equation solved in 
Ref. [2(| is a time- local differential equation and it can be 
obtained from cq. ([2U)) by replacing p(t') with p(t). 

IV. COMPARISON BETWEEN EXACT AND 
APPROXIMATE SOLUTIONS 

Exploiting the results obtained in the previous sections 
we now analyze the time behavior of the two qubits com- 
paring in particular the exact dynamics with the ones 
stemming from the NZ and TCL approaches. Our aim 
is to highlight the performances of two approximate ap- 
proaches and to point out their range of validity. As un- 
derlined in the introduction there is not indeed a general 
theory that predict which one of the two methods is to be 
preferred and generally speaking their range of validity 
strongly depends on the specific features of the system, 
namely the interaction hamiltonian, the interaction time, 
the environmental state and the spectral density. 

In what follows we assume that each qubit interacts 
resonantly with a reservoir with Lorentzian spectral den- 
sity 

Mu) = J 3 (w) = J(v) = ±- , (27) 

2ir [ujo — + y 

where 7 is a parameter which in the Markovian limit co- 
incides with the system decay rate, and A is the reservoir 
bandwidth. This is for instance the case of two atoms 
interacting each of them with their own cavity field in 
presence of cavity losses Thus the two correlation 

functions f\(t — t') and fi{t — t') defined by eq. (fl"5|) co- 
incide and are given by 

hit = f 2 (t = f(t = | T Ae- A l«-*'l. (28) 




FIG. 1: Time evolution of the population of the state 1 10) for 
a system initially prepared in the Bell state (|10) — |01)). 
The width of the Lorentzian spectral density is A = IO7, the 
strength of the coupling constant between the two qubits is 

n = 0.0017. 



From the latter equation, it is clear that the bandwidth 
A plays the role of the inverse of the reservoir memory 
time. 

The system dynamics will be analyzed considering the 
two baths in a thermal state at T = whereas the 
two qubits are supposed at t = in the Bell state 
\4>i0)) = -j= (| 10) — 1 01 ) ) . Let us concentrate on the pop- 
ulation Pio(i), that is on the probability of finding the 
two qubits in the state |10). To test the validity of the ap- 
proximate approaches we explore three different regimes 
varying the width of the Lorentzian spectral density A. 
This investigation will allow us to assess in which cases 
the solutions of the master equations are efficient in the 
description of the true dynamics of the system. 

Figure |T]) shows a comparison among the exact, the 
TCL and the NZ solutions in the case of large reservoir 
bandwidth A = IO7. The plot is done against the dimcn- 
sionless variable jt and the coupling constant between 
the qubits is fixed to f2 = O.OOI7. We can clearly appre- 
ciate the perfect agreement between the exact analytical 
solution and the approximate ones for the short time be- 
havior but also for long interaction times. In this case the 
two approaches TCL and NZ both provide a very good 
description of the dynamics and we may conclude that 
there is no way to establish if one method is to be prefer- 
able with respect to the other, they indeed give the same 
results. However in such cases the TCL master equation 
might be preferred since it involves a time local first or- 
der differential equation and therefore it is easier to solve. 
In the inset of fig. ((Tj) the short time behavior is shown. 
It is interesting to underline the initially quadratic be- 
havior that witnesses the non-Markovian features of the 
dynamics. 

In fig.© the same quantity is reported choosing A = 7. 
Despite the good agreement for the short time dynam- 
ics, we observe significant deviations when time increases. 
In particular concerning the long time behavior, the NZ 
equation leads to a very bad approximation. For times 
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the exact evolution of the populations Pio (t) are not due 
to the spin-spin interaction constant f2, which we have 
taken small on purpose. The oscillations are instead due 
to the fact that, in this case, the Lorentzian peak is so 
narrow to make the environment equivalent to a cavity 
with losses, as one might verify for example by using the 
pseudomode approach in Ref s . [TTI. H3, [HI. Il8j . Therefore, 
our exact solution in this regime, could be exploitable to 
describe the dynamics of two interacting qubits put in- 
side two different optical cavities. 



FIG. 2: Time evolution of the population of the state 1 10) for 
a system initially prepared in the Bell state (|10) — |01}). 
The width of the Lorentzian spectral density is A = 7, the 
strength of the coupling constant between the two qubits is 

n = 0.0017. 
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FIG. 3: Time evolution of the population of the state |10) for 
a system initially prepared in the Bell state (|10) — [01)) . 
The width of the Lorentzian spectral density is A = O.OI7, 
the strength of the coupling constant between the two qubits 

is n = 0.0017. 



longer than some critical values the solution for the popu- 
lation Piq (t) cannot represent a true diagonal element of 
a density matrix anymore, since it indeed assumes neg- 
ative values. We may conclude that for this range of 
parameters the TCL solution gives a better description 
of the dynamics since it reproduces all the qualitative 
features of the exact solution. 

Finally in fig. ((3]) we examine the regime A = 0.0 17 
which, according to eq. (P5)) . corresponds to very strong 
reservoir correlations and very long memory time. We 
observe once again a perfect agreement among all the 
three approaches in the short time behavior but in this 
case we assist at a failure of the TCL approach too. 
The solution of the TCL master equation (dashed line) 
doesn't succeed to follow the Rabi's oscillations witnessed 
by the exact dynamics (solid line). The NZ approaches 
presents the same problem of not conserving the positiv- 
ity of the density matrix as in the previous case. Thus, 
in this case, both the perturbativc approaches are not 
suitable to describe the dynamics of the system. It is 
interesting to observe that the oscillations appearing in 



V. DISCUSSION AND CONCLUSIVE 
REMARKS 

The problem of the proper description of open quan- 
tum systems is still far from having a complete and gen- 
eral solution. In particular, while the features of the 
Markovian dissipative dynamics are all well established 
and accepted, a lot of work has still to be done to claim 
general statements on the validity of the different possible 
non-Markovian approaches available in the literature. 

In this paper we concentrated on the study of two pop- 
ular methods aimed at describing non-Markovian dynam- 
ics: the Nakajima-Zwanzig and the time-convolutionless 
master equation approaches. Exploiting an exact solu- 
tion for the dissipative dynamics of two coupled qubits 
interacting with independent reservoirs, we have shown 
that the TCL approach reproduces all the features of the 
non-Markovian dynamics for a range of parameters much 
wider than the one in which the NZ equation gives results 
which are physically reasonable, since the latter approach 
may violate the positivity condition on the density ma- 
trix already for reservoir correlations which are not very 
strong. 

The discrimination of the best master equation ap- 
proach for the problem under study is a very important 
issue, because the problem of the dissipative dynamics 
of two interacting qubits has been given a lot of atten- 
tion in recent years, especially from the point of view 
of the entanglement dynamics [H, GJ, HI US US [M|. 
So, once we have shown that for some parameters the 
TCL master equations provides a very good description 
of the non-Markovian dynamics, one may use the same 
approach for the description of temperature effects on the 
system dynamics, since we do not have at our disposal 
an exact solution of the dynamics at general reservoir 
temperatures. 

We finally note that our exact solution can moreover 
be exploited in order to study the zero temperature en- 
tanglement dynamics in a range of parameters in which 
none of the two master equation approaches can be used: 
in this way one could extend the studies given in Refs. 
[3 LIS, Llo| , so that one could get the most general fea- 
tures of the quantum correlations between two qubits 
coupled to two independent reservoirs. These points will 
be the subject of our future research. 
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